Decohering d-dimensional quantum resistance 

Dibyendu Roy, N. KumaiQ 
Raman Research Institute, Bangalore 560080, India 
(Dated: February 1, 2008) 

Abstract 

The Landauer scattering approach to 4-probe resistance is revisited for the case of a d- 
dimensional disordered resistor in the presence of decoherence. Our treatment is based on an 
invariant-embedding equation for the evolution of the coherent reflection amplitude coefficient in 
the length of a 1-dimensional disordered conductor, where decoherence is introduced at par with 
the disorder through an outcoupling, or stochastic absorption, of the wave amplitude into side 
(transverse) channels, and its subsequent incoherent re-injection into the conductor. This is essen- 
tially in the spirit of Biittiker's reservoir-induced decoherence. The resulting evolution equation for 
the probability density of the 4-probe resistance in the presence of decoherence is then generalised 
from the 1-dimensional to the d-dimensional case following an anisotropic Migdal-Kadanoff-type 
procedure and analysed. The anisotropy, namely that the disorder evolves in one arbitrarily chosen 
direction only, is the main approximation here that makes the analytical treatment possible. A 
qualitatively new result is that arbitrarily small decoherence reduces the localisation-delocalisation 
transition to a crossover making resistance moments of all orders finite. 
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I. INTRODUCTION 

Electron localisation , strong as well as weak, and the associated metal-insulator 
transition and conductance fluctuations [2I are due essentially to the time-persistent in- 
terference of the complex wave amplitudes that result from multiple elastic scattering on 
randomly distributed defects in the conductor with quenched potential disorder. Similar in- 
terference effects also manifest in mesoscopic systems as various phase-sensitive phenomena, 
e.g. the well known persistent ring currents and the Aharonov-Bohm oscillations. Clearly, 
these one-electron phase-sensitive phenomena can get suppressed by decoherence. Micro- 
scopically 4|, decoherence can arise from incoherent processes involving, e.g., the inelastic 
electron-phonon or the electron-electron scattering, as also from an entanglement with the 
environmental degrees of freedom that remain undetected or unmeasured. (While coher- 
ent inelastic scattering is in principle possible, as indeed in the case of neutron scattering, 
it is not relevant to the case of coherent multiple scattering of electrons in a disordered 
conductor). The question now is how to incorporate decoherence phenomenologically in 
an analytical treatment of the otherwise Hamiltonian system such as the system of non- 
interacting electrons moving in a lattice, or a continuum with random elastic scatterers, 
e.g., the Anderson model system for metal-insulator transition in random lattices [lj. De- 
coherence has often been included theoretically and probed experimentally through a phase 
breaking, or dephasing cut-off length scale introduced on physical grounds {2], {J. It is clearly 
desirable, however, to have a phenomenology for introducing the degree of decoherence in 
the analytical treatment of elastic scattering in a disordered conductor. A highly successful 
and widely used approach to decoherence was pioneered by Biittiker through 
the idea of reservoir-induced decoherence. The latter could be introduced naturally in the 
scattering approach of Landauer [lOj to quantum transport, e.g., the 4-probe resistance. 
For the reservoir- induced decoherence, one inserts a scattering matrix with appropriately 
chosen side (transverse) channels, and thereby outcouple a partial wave amplitude into an 
electron reservoir. The amplitude re-emitted from the reservoir is then re-injected back into 
the conductor, adding necessarily incoherently to the transmitted amplitude along the con- 
ductor (the longitudinal channel) that carries the transport current. The chemical potential 
of the reservoir is, of course, tuned so as to make the net current in the side channel vanish 
on the average. (This is clearly analogous to the "potentiometric" probe of Engquist and 
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Anderson The net result is the introduction of decoherence, or partial coherence, that 

can be readily parametrised. It describes, in particular, the quantum-to-classical crossover 
of a series combination of conductors 7] with increasing strength of the coupling to the in- 
tervenin g re servoirs. While used extensively in the context of mesoscopic (zero-dimensional) 
systems [121 ]. the reservoir- induced decoherence has also been invoked by many workers for 
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partial coherence in quantum transport on tight-binding lattices - without disorder 



151 ]. and with weak disorder 
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Il8j |. as also in a disordered continuum |l9 |. 
These studies are, however, confined to 1-dimensional conductors. 

In this work, we have considered the case of a d-dimensional conductor for d > 1 in the 
presence of both the quenched disorder and decoherence. Our analytical treatment is based 
on the invariant-embeddin g ap proach developed earlier for a 1-dimensional conductor with 



quenched disorder [20j, 
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221 ]. and its subsequent generalisation to higher dimensions for 
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24j . Here, first the elastic scat- 



anisotropic disorder using the Migdal-Kadanoff technique 
terers (resistances) are combined in series quantum-mechanically along an arbitrarily chosen 
direction and then classical Ohm's law is used to combine these resistances in parallel along 
the transverse directions. This is followed by a scaling transformation with an infinitesimal 
increase in scale at each step. The resulting 'transverse' mixing up of disorder through the 
evolution equation is known to give a qualitatively correct description of the weak scattering 
regime in the absence of decoherence [24|, despite the assumption of anisotropic disorder, 
which is an approximation. In our approach, decoherence and elastic scattering (quenched 
disorder) are treated formally at par through a proper insertion of the scattering (S-) ma- 
trices, i.e., transverse channels distributed over the conductor. Specifically, a side-channel 
is to be viewed as causing a stochastic absorption - a coherent process by itself. The in- 
coherent re-injection with zero net side-current is, however, effectively realised through the 
use of the Landauer expression |i?(L)| 2 /(l — \R(L)\ 2 ) for the 4-probe resistance, but with 
|i?(L)| 2 now calculated as the coherent-only reflection coefficient. A physically robust argu- 
ment is presented for the self-consistency of this procedure. The main results derived are, 
(a) elimination of the metal-insulator transition (the unstable fixed point) for an arbitrarily 
small strength of decoherence. This is indeed expected on physical grounds inasmuch as 
metal-insulator transition with increasing static disorder is essentially due to the coherent- 
back scattering [2] (where the back scattered amplitudes traversing the time-reversed paths 
add up in phase), while decoherence suppresses this phase coherent effect; (b) suppression 
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of the 4-probe resistance fluctuations with increasing decoherence strength making all the 
resistance moments finite; and (c) a correction to conductivity due to decoherence in the 
metallic limit that mimics the conventional phase cut-off length scale. 



II. MODEL AND INVARIANT-EMBEDDING: 1-DIMENSIONAL CASE 

Consider a model Hamiltonian H for the system of non-interacting electrons in a 1- 
dimensional disordered conductor of length L: 

ft 2 d 2 

H = -^^- 2 +V(x), (2.1) 
2m ox z 

where V(x), < x < L is a spatially random potential (quenched disorder) assumed to be 
delta-correlated Gaussian as 

(V(x)V(x)} = V 2 5(x-x') . 

Let an electron wave of unit amplitude be incident at Fermi energy (E F = h 2 k 2 F /2m) on the 
sample from right, and let R(L) and T(L), respectively, be the reflection and the transmission 
amplitude coefficients. Next, let the sample of length L be embedded invariantly in a 
supersample of length L + AL (Fig. [TJ. It is readily seen that the elastic scattering from 
the random potential in the interval AL with kpAL << 1 can be viewed as due to a delta- 
function potential of strength V(L)AL, the corresponding scattering matrix being ASe 

( 2mVAL i , 2mVAL \ 
2itfk F 2ifPk F 
i i 2mVAL 2mVAL J v ' 

1 2ih 2 k F 2ih?k F ) 

This gives an evolution equation for the S'-matrix in the sample length L. Specifically, we 
have for the amplitude reflection coefficient 
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^ = z^(L)(l + R(L)) 2 + 2ik F R{L) , (2.3) 
with £(L) = and (^(L)^(L')) = A5(L — L'). 

We are now in a position to introduce decoherence at par with the random elastic scattering 
within this approach. We recall the 4x4 S'-matrix with the side channels as introduced 
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FIG. 1: Shows disordered sample of length L embedded invariantly in a supersample of length 
L + AL. Shown also are the elementary matrices for the elastic (ASe) and the decohering (ASd) 
scatterings in AL, with the incident, the transmitted, and the reflected waves at Fermi wavevector 

by Biittiker Q: 
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Here the outcoupling through the side channels is parametrised by e, which must be of order 
AL in the present case. Accordingly, we use the 2x2 sub-matrix 





AS 



D 




(2.5) 



for insertion into the interval AL. It describes the outcoupling into the side channels, i.e., 
the stochastic absorption, as also the coherent transmission directly through the interval AL. 
(Its connection with the reservoir-induced decoherence will be clarified below later). Figure 
(CQ) is a schematic depicting the insertion of the elementary ASe and ASd in the interval 
AL. Clearly, for A^AL << 1, the exact spatial order and the locations of the two insertions 
within the interval AL are not relevant. Combining these two elementary S-matrices (ASe 
and ASd) for AL with the S— matrix (S(L)) for the sample of length L in series, we can 
read off the emergent quantities R(L) and T(L): 

AT 2 e 2lkFAL R(L) 



R(L + AL) = AR + 



1 - AR R(L) e 2ik F AL 



(2.6) 



with AR 



2mVAL 
2ih 2 k F 



and AT 



e + 



2m^AL 
ih 2 kF 



In the limit AL — > 0, we obtain the evolution equations for the amplitude reflec- 
tion/transmission coefficients R(L) and T(L): 

^ = z^(L)(l + R C (L)) 2 + 2ik F R c {L) - V R C (L) , (2.7) 
and ^ = iy£(L)(l + Rc(L))T c (L) + ik F T c (L) - ^T C (L) , (2.8) 

where rj = e/AL, AL — > parametrises decoherence. Here we have introduced the 
subscript 'c' just to emphasize that the reflection/transmission amplitude coefficients in 
Eg. (12.71) are coherent. 

It seems in order at this stage to clarify how decoherence is realised in relation to the 
sample resistance by the insertion of the side channel through ASp. Clearly, the embedding 
Eq. (12.71) describes evolution of the coherent reflection amplitude R C (L). (Similarly, T C (L) 
is the coherent transmission amplitude as depicted in Fig. ([I]). The embedding equation 
for T C (L), however, is not autonomous - it involves R JL )). The outcoupling into the 



side channels corresponds to a stochastic absorption [25|, [26|, l27j in the interval AL. This, 
however, has to be re-injected now incoherently back into the conductor. Inasmuch as this re- 
injected current necessarily flows down the chemical potential gradient, it contributes to the 
total transmitted current equal to (within constant of proportionality) \T C (L)\ 2 + \T in (L)\ 2 = 
\T tot (L)\ 2 , where the subscript 'in' denotes incoherent. From the conservation of the total 
current flowing down the conductor, we must have |T C (L)| 2 + \T in (L)\ 2 = 1 — |_R C (L)| 2 . 
Now, recall that the Landauer resistance (p^' D ^) formula p( d ' D ) = (1 — \T tot \ 2 ) /\T to t\ 2 holds 
for arbitrary |T tot | 2 (coherent or incoherent both). Here the superscript {d,D) denotes the 
dimensionality d and the decoherence parameter D. Thus, we have p^ d ' D ^> = |i? c | 2 /(l — |i? c | 2 ) 
given entirely in term of R C {L) which is calculable from Eq. (12.71) . Thus the 4-probe resistance 
|.R C (L)| 2 /(1 — |i? c (L)| 2 ) incorporates self-consistently the incoherent re-injection. Here we 
must re-emphasize that \R C (L)\ 2 is the coherent reflection coefficient given by and calculable 
from the embedding Eq.( l2.7j) . 

Our next step is to obtain the 'Fokker-Planck' equation for the probability density of the 
reflection coefficient r(L) = \R C (L)\ 2 from the stochastic differential E q. (12.71) w hich serves 



as the Langevin equation here. Following the now familiar procedure 
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28|], we 
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obtain 

= ^^(^^(D^^+^^pd)^^ , (2 . 9) 

L 2 
with I = — l — ——5- and D = 2r?Z . 

This is clearly a two-parameter (7 and D) evolution equation. The two independent pa- 
rameters Iq and D are, of course, composed of the two basic independent parameters Akp 
(measure of disorder) and rj (measure of decoherence) . Thus, e.g., D may vary through 77, 
while Iq can remain constant. 

Equation (12.91) in the limit of large length L » l gives a steady-state distribution Poo(r) 
for the reflection coefficient r 



p M (r) = |£|exp(|D|)exp(-f 
fl — r " 



n_,v2 ' '• < L ( 2 - 10 ) 



Note that for D = 0, the limiting distribution in Eq. (12.101) tends to the delta- function, 5(1 — 
r), and not to zero. (This can be readily seen by noting that the probability distribution is 
normalized to unity for all D). This means that the reflection coefficient becomes unity with 
probability one, as it must for an infinitely long Id disordered wire without decoherence (all 
states being localized then, a well known result from Anderson localization in one dimension). 
The corresponding resistance moments are all finite for D ^ 0. In particular, the limiting 
value of the average 4-probe resistance in the presence of decoherence is 



nh r 



e z 1 — r 



e 2 \D\ ' 

With this preparation (Eq. (12.9[) ) in hand, we now turn to the case of d-dimensions. 



(2.11) 



III. HIGHER-DIMENSIONAL CASE 

Changing over to the 4-probe resistance p = r/(l — r) (measured in the unit of 7rh/e 2 ) as 
the new variable with the associated probability density P^(p,l), Eq. (I2.9p reduces to 

-^T = p(p + l )~^r + {( 2 p + !) + d p(p + 1)}^- + d ^p + l )P {1) ■ (3-i) 

The corresponding nth resistance moment in 1 dimension is 



rn 



poo 

/ pM(p,l)p n dp . (3.2) 
Jo 



Multiplying both sides of Eq. fl3.ip by p n and integrating by parts on the RHS, we get the 
evolution equation for the 1- dimensional moment 

^ = n(n + l)pS^ + nVK } - Dnp^ - Dnpffi , (3.3) 

which is hierarchical in nature (i.e., the equation for p«' D ^ involves Pn^P and Pn+P)- For 
D — 0, however, the equation for pn involves the lower-order moments only leading to 
a closure of the hierarchy. Thus, the presence of decoherence (D ^ 0) brings about a 
qualitative change in the structure of the coupled equations for the moments of different 
orders. For D = 0, the solutions of Eq. (l3.3p for the 1st and the 2nd moments are readily 
obtained as 

p\ m = \^ - 1) , 

In writing the last equation above, we have eliminated the length I in favour of an implicit 
relation between pj, 1 ' '' an d Pi 1 ' ' '■ We have verified by iteration of Eq. (13.31) . that this relation 
remains valid for p^' D ^ and p[ 1 ' D ' ) to a good approximation for D ^ 0, and will be used as 
such. Substituting for p^' D ^ in terms of Pi' D ^ in Eq. (13. 3p for n = 1, and integrating we 
obtain a relation between I and p[ 1 ' D ' ) 

1 = Jo -lDpf' D)S - 2Dp i l' D)2 + (2 - D)p^ D) + 1 ' (3 ' 5) 
(From now dummy integration variable will be distinguished by a tilde). Hereinafter, the 
superscript D in Pi' D ^ will be dropped except when required for the sake of clarity. Defining 
the associated moment generating function \ (^jO an d the cumulant generating function 
K<U(x,l) of P (1) (p,/) as 

POO 

X W (x,l) = / e- x <>pW{p,l)dp , 
Jo 

KW(x,l) = \n X {1) (x,l) , 

we derive from Eq. (l3.3p their evolution equations 

\- = (x 2 + Dx)^_ + (2x -Dx-x 2 )^- , (3.6) 

—— = (x 2 + Dx) 09 + (x 2 + Dx)(^—) 2 + (2x - Dx - x 2 ) — x. 3.7 

at ox 2 ox ox 



Now, we proceed to generalise the above equations to the case d > 1 . For this we closely fol- 



low the Migdal-Kadanoff procedure as in Ref. 



23] , assuming the quenched disorder to evolve 



along one chosen direction only. This anisotropic disorder is admittedly an approximation, 
but it is known to reproduce correctly the qualitative features of the Anderson transition in 
the absence of decoherence, as shown in the earlier works 23 



241 ] . The probability density 
(p, I) of the resistance of a d-dimensional hypercubic sample is accordingly found to obey 
the integro-differential evolution equations 



d x (d) 

dlnl 
dK^ 
dlnl 



-{d - l)x— \ 



pp. (d) 

Dx) — \- + {2x-Dx 



dx 



(x 2 + Dx) 



dx 2 
d 2 K^ 
dx 2 



x 



Dx){ 



d X {d) 
dx 
dK^ 2 
dx 



x X {d) ] I , 

(2x — Dx — x 2 ] 



x) I , 



(3.8) 

dx 
(3.9) 



where / in the above equations is given by the integral in Eq. fl3.5p . but with in the 
integrand now re-interpreted as pf 1 . Clearly, in the limit D — 0, the above equations for 
the generating functions reduce to the corresponding Eqs.(6,7) of Ref. 2^ |. 

In particular the fixed point probability distribution for d — 3 obtained by setting 
dx^/dhi I = and inverting the Laplace transform of the solution for nothing 
but the known fixed point power-law distribution 24]. 

In the presence of decoherence (fl / 0), however, there is no fixed point even for arbi- 
trarily small values of D for d = 3. In order to see this, consider the evolution equation for 
the first cumulant k[ (= p±) obtained from the cumulant-generating Eq. (13.91) 

dK[ d) 



dlnl 



(d - l)K[ d) + [1 + 2K {d) - 
dp? 



DK {d) - DK\ 



(dy 



DK. 



(«0i 



(3.10) 



-ppf Y - 2Dpf )2 + (2 - D)pf ) + 1 ' 

where we have replaced the length I in terms of p^ as explained above. Carrying out the 
integration occurring in Eq. fl3.10l) numerically (using Mathematica) , we found no solution 
with dK? /dlnl = for any non-zero value of D however small (down to D ~ 10~ 6 ) 
confirming that there is no fixed point. This should, of course, be physically so inasmuch 
as the decoherence is expected to suppress quantum interference effects (and localisation), 
in the limit of large sample size. For D ^ 0, however, we do expect the probability density 
to vary slowly in the vicinity of the D = fixed point, now become a crossover. Indeed, 
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setting dx^ /dial ~ for small non-zero D, we obtain for the quasi-fixed-point probability 
density of resistance 

P(pf) = ~ [ J Pl } , (3-11) 



T(l -a,D) 



where r(l — a, D) = e u u a du and a 



D t\ p W 



Here p± (~ p± = 1.96 for d — 3) is the average resistance corresponding to the quasi- 
fixed-point probability density, and l\ <d) is the value of the integral Eq. fl3.5p with the upper 
limit pf\ It is clear from Eq. fl3.lll) that a non-zero value of D (decoherence) makes all the 
resistance moments finite, that is it cuts-off the otherwise divergent resistance fluctuations. 
In the absence of decoherence (D=0), Eq. (13.1 11) gives a power -law probability distribution 
for resistance at the mobility edge as in Ref . [24| . It is to be noted, however, that all numerical 
work on tight-binding Anderson model shows that the distribution of conductance in 3d at 
the mobility edge (the fixed point) is far from a power law 29_|, |30(]. We think that this 
could be for two reasons: First, the neglect of the transverse fluctuations in our anisotropic 
Migdal-Kadanoff procedure, and secondly as the numerical results are all for the ensemble 
averaged two-probe conductance (Tr tP) (where t is the transmission matrix) while we have 
calculated the ensemble averaged four-probe resistance. It is to be noted here that while the 
four-probe resistance is unbounded from above and can, therefore, have large fluctuations, 
the two-probe conductance is by definition bounded from above and can fluctuate relatively 
much less. Thus, e.g., in the Id case, Tr tv < 1; but, of course, there is no fixed point in the 
Id case. We would like to point out here that the invariant imbedding equation is, of course, 
known for d-dimensional as also for the quasi-one dimnsional case (see Ref. 22], Eq.2.28), 
but an analytical solution is lacking, making any comparison with the available results for 



the quasi- Id, D=0 conductance distribution [3l|, |32J impossible. Our immediate interest, 



however, lies in the fact that decoherence cuts off the resistance fluctuations exponentially. 

Finally, we consider the asymptotic behaviour of the resistance in 3 dimensions in the 
presence of decoherence in the metallic regime as the sample size tends to infinity. In 3 
dimensions with D ^ 0, we expect the resistance to tend to a small value in the mean along 
with a narrow width (the variance) of the distribution. This motivates us to approximate 
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the evolution Eg. (13. 10p for the first moment as 
dp? 



(d-l)p? + [l + (2- D)pl 



dlnl 



(d) 
Pi 



dpi 



(3.12) 



/o 1 + (2 - D)pf ' 

Now, consider first the 3-dimensional case (d=3) in the metallic regime starting with the 
resistance p = pf\lo) at a length scale l . Let this evolve through Eq. fl3.12l) to a length 
scale I >> Iq with p\ (I) = p « p . Eq. fl3.12l) then gives 

dp? 



p[ 3) 



~(3) 2-D ~(3) 2 
po ~Pi H —Pi 



ln( 



(3.13) 



or, in term of the conductivity a^il) = g/l, g = 1/p and go = 1/po, we have 

r (3)/,A 



(3.14) 



Equation (13.141) clearly shows that increasing decoherence (D) increases the metallic con- 
ductivity in 3 dimensions. Indeed, one can re- write the correction D/21q as l/L^ with 
a phase-cut-off (dephasing) length scale as usual. Proceeding in similar way, we get for 
the 2-dimensional case a logarithmic correction to the conductivity (I) (noting that in 2 
dimensions conductivity is the same as conductance) 



^ (2) (0 = ^ + ^in(f; 



(3.15) 



where <To is the conductivity (or the conductance) at the starting length scale Iq. Again, the 
conductivity cr^il) is seen to increase with increasing decoherence D. This is qualitatively 
consistent with the negative temperature coefficient of resistance observed in disordered 
conductors at low temperatures in the weak localisation regime, in particular for 2d systems 

0. 

IV. DISCUSSION 

We have extended the phenomenology of decoherence known well in the context of 
)hase-sensative systems, such as mesoscopic rings [33j and 1-dimensional quantum wires 



3. 
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m 



13, 



181 ]. to higher dimensions - specifically to a d-dimensional disordered 



conductor for d = 2 and 3. Our treatment here follows the invariant embedding approach 
developed earlier l2ll . |22|, beginning with the Id case. It treats decoherence and dis- 



order formally at par in that the two are introduced through appropriately chosen and 
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parametrised scattering matrices distributed over the conductor. The invariant embedding 
appraoch is naturally suited to the problem on hand as it gives the evolution-in-length 
of the resultant emergent quantities such as the reflection coefficient related directly to 
the Landauer 4-probe resistance of interest. Decoherence is realised specifically through 
stochastic absorption of the wave-amplitude into distributed side (transverse) channels, and 
the subsequent re-injection of the absorbed fraction back into the conductor so as to add 
incoherently to the (longitudinal) coherent transport. This is essentialy in the spirit of 
Biittiker's reservoir-induced-decoherence. A point to note here is that the current-conserving 
re-injection is realised here self-consistently through the use of the 4-probe resistance which 
now needs to be calculated with the coherent-only reflection coefficient. Extension to higher 
dimensions has been carried out within the Migdal-Kadanoff procedure assuming the dis- 
order to evolve only along an arbitrarily chosen direction for the current. This choice of 
anisotropic disorder is admittedly an approximation, but its innocuous nature is borne out 
a posteriori by the fact that this approximation had correctly given the unstable fixed point 
for the disorder induced Anderson (metal-insulator) transition for d = 3 in the absence of 
decoherence. Its reasonableness may be attributed to the transverse mixing up of disorder 
by the evolution equation. Physically, moreover, the classicalisation expected from decoher- 
ence should make the approximation even better. A non-trivial result of our work is the 
elimination of the unstable (Anderson) fixed point due to decoherence. Again, it is expected 
on physical grounds that the fixed point should get replaced by a crossover for D ^ 0. So 
is the finiteness of all moments, that is the suppression of resistance fluctuations due to 
decoherence, as is evident from our Eq.f l3.lll) . A point to note is the decoherence correction 
to the quantum conductivity for d = 3, where a cut-off length (dephasing length) appears 
naturally. Finally, we would like to point out here that the decoherence, through stochas- 
tic absorption into the transverse channels and the re-injection, does not cause scattering 
in the coherent longitudinal (transport) channel in the sense of momentum randomisation 
that would have given additional resistance. Indeed, as is clear from our Eq. fl2.7p . in the 
absence of scattering by disorder, the reflection amplitude (R) remains identically zero for 
all lengths, independently of the value of rj (that parametrises decoherence). This is also 
obvious from the Eq. (12.51) . We would aptly like to call this a pure decoherence without any 
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concomitant elastic scattering. 
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